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AN  INTEGRAL  EQUATION  FOR 
THE  FLOATING  BODY  PROBLEM 


by 

T.  S.  Angel 1,  G.  C.  Hsiao,  R.  E.  Kleinman 
Mathematical  Sciences 
University  of  Delaware 
Newark,  Delaware  19716 


y  The  time  harmonic  three  dimensional  finite  depth 
floating  body  problem  is  reformulated  as  a  boundary 
integral  equation.  Using  the  elementary  fundamental 
solution  that  satisfies  the  boundary  condition  on  the 
sea  bottom  but  not  the  linearized  free  surface  condition, 
the  integral  equation  extends  over  both  the  ship  hull 
and  the  free  surface.  It  is  shown  that  this  integral 
equation  is  free  of  irregular  frequencies,  that  is,  it 

has  at  most  one  solution.  - — 

1.  Introduction 

In  his  classic  work  on  the  floating  body  problem,  F. 

John  (1950) ,  showed  how  the  boundary  value  problem 
could  be  reduced  to  an  integral  equation  over  the  wetted 
portion  of  the  ship  hull.  The  kernel  of  his  integral 
operator  was  the  Green's  function  for  the  entire  fluid 
domain  with  no  ship  present  that  satisfied  the  boundary 
condition  at  the  bottom  of  fluid  (assumed  flat)  and  the 
linearized  free  surface  condition  on  the  entire  fluid-air 


r.  •' 


boundary.  John  demonstrated  the  existence  of  irregular 
frequencies ,  frequencies  for  which  the  integral  equation 
was  not  uniquely  solvable.  Recently  Kleinman  (1982) 
provided  two  methods  of  modifying  the  integral  equation 
so  that  there  were  no  irregular  frequencies.  In  one 
case  the  domain  of  the  integral  operator  was  enlarged  and 
in  the  other  the  operator  itself  changed,  but  both  methods 
employed  John's  Green's  function  which  is  rather  complicated 
especially  in  the  three  dimensional,  finite  depth  case. 

Another  way  to  treat  this  problem  is  to  employ  a  much 
simpler  Green's  function,  one  that  satisfies  only  the 
boundary  condition  at  the  bottom  of  the  fluid.  Since  this 
does  not  satisfy  the  free  surface  condition,  we  get  an 
integral  equation  over  both  the  wetted  surface  of  the  ship 
hull  and  the  free  surface.  Such  an  integral  equation  has 
been  derived  and  even  solved  numerically  for  certain  cases, 
e.g.  Yeung  (1978)  and  Bai  and  Yeung  (1974) .  Numerical 
evidence  indicated  that  this  integral  equation  did  not 
exhibit  irregular  frequencies  but  no  conclusive  analytical 
argument  has  yet  appeared  to  support  this  conjecture. 

The  present  paper  provides  a  proof  of  the  conjecture 
that  this  integral  equation  has  no  irregular  frequencies. 


By  irregular  frequencies  is  meant  frequencies  for  which 
the  integral  equation  is  not  uniquely  solvable  even 
though  the  solution  of  the  corresponding  boundary  value 
problem  is  unique.  What  we  prove  is  that  the  integral 
equation  obtained  using  a  simple  combination  of  elementary 
sources  is  uniquely  solvable  at  all  frequencies. 

It  should  be  emphasized  that  our  concern  here  is 
not  with  uniqueness  for  the  boundary  value  problem  itself. 
There  John  required  certain  geometric  restrictions  in 
order  to  establish  uniqueness.  These  may  be  somewhat 
relaxed  to  include  hull  forms  with  comers  and  non  normal 
intersections  with  the  free  surface  (see  Kleinman,  1982) . 
However,  in  the  three  dimensional  case  treated  here,  we 
retain  the  restriction  that  vertical  rays  from  the  free 
surface  may  not  intersect  the  ship  hull  in  order  that  the 
boundary  value  problem  be  uniquely  solvable.  Our  concern 
here  is  with  integral  equation  formulations  and  the  irregular 
frequencies  which  are  introduced  in  some  cases. 

It  should  be  noted  that  the  occurrence  of  irregular 
frequencies  in  integral  equation  formulations  of  acoustic 
scattering  problems  is  entirely  analogous  to  the  present 
case.  (See  e.g.  Smirnov,  1964;  Brundrit,  1965;  Copley,  1968; 
Schenck,  1968  and  Chertok,  1970,  1971)  However,  methods 
for  removing  the  irregular  frequencies  in  acoustic  scattering 
all  essentially  involve  making  the  kernel  of  the  integral 
equation  more  complicated  (e.g.  Brakhage  and  Werner,  1965, 
Burton  and  Miller,  1971,  Kleinman  and  Roach,  1974,  1982) . 
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Specifically,  we  treat  the  three  dimensional  floating 
body  problem  with  finite  depth,  h.  If  we  denote  the  fluid 
domain  by  D+,  the  hull  by  CQ,  the  free  surface  by  and 
the  bottom  by  Cfi,  and  if  we  denote  by  D_  the  domain  consisting 
of  the  upper  half  space  and  the  interior  of  the  ship  hull, 
then  geometry  may  be  illustrated  as  in  Figure  1. 


Figure  1 

The  function  $  solves  the  floating  body  problem  if 

V2,  .  o  in  D+,  |±  -  V  on  Cc,  |i  -  0  on  CB  ,  (1 

.|n  +  k$  *  0  on  Cf  • 

and  provided  $  satisfies  a  radiation  condition.  Here 

an 

is  the  normal  derivative  directed  into  D+  and  V  is  a  given 
function.  The  radiation  condition  is  specified  in  the  form 


44-  -  ik^$  ■  o  (p-1/2)  as 


P  -►  CO 


(2 
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uniformly  in  6  and  y.  This  condition  may  be  shown  to 
guarantee  that 

eikop  -1 

♦(p)  ■  - — —  (f(e)  +  o<p  x) )  as  p  +  •  ,  (3) 

Sp 

(p , 8 )  being  polar  coordinates  in  the  free  surface-water 
plane  and  kQ  is  the  positive  real  of  the  transcendental 
equation 

k  =  k  tanh  k  h  .  (4) 

o  o 

Now  define  the  Green's  function 


Y  (p#q) 


1  1 
MP=qT  ”  2tt  |p-q1 


(5) 


where  p 


(xp»  yp>  2p>'  <3 


(Xq,  -2h-yq, 


(Xq»  yq.  *q)  and  qx 
and  we  have  oriented  a  rectangular  coordinate  system  so  that 
the  plane  y  -  0  is  the  water  plane  and  free  surface  while 
y  -  -h  is  the  bottom. 

With  the  Green's  function  defined  in  (5),  which 
has  a  double  strength  singularity  on  Green's  theorem 
for  solutions  of  Laplace's  equation  in  D+  which  satisfies 
the  radiation  condition  (2)  takes  the  form 


CouCf uCB 


(yCp.  q)  !“ - $<q)  )  ds  *  o(p)  4>  (p) 

q  q  H 


(6) 


where 


a(p)  ■  2  for  p  €  D  u  C 

»  1  for  p  on  smooth  points  of  CQ  u  Cj  (7) 

■  0  for  p  £  D_. 

If  ♦  satisfies  all  of  the  boundary  conditions  in  (1)  we  obtain  the 
boundary  integral  equation 


*  /  Y(p,q)V(q)ds 
r  “ 


(8) 


where  p  lies  either  on  CQ  or  C^.  The  integral  on  CB  vanishes 

since  both  y  and  <p  satisfy  a  homogeneous  Neumann  condition 

and  the  integral  over  a  large  cylinder  vanishes  since  y  » 

-1  -1/2 

0 (P  )  and  <J>  ■  0(p  ),  the  radiation  condition  ensuring 

that  <J>  has  asymptotic  growth  given  by  (3).  As  explained  in  the 
introduction,  this  equation  has  irregular  frequencies  if  there  are 
certain  values  of  k  for  which  the  homogeneous  equation  (V=0) 
has  nontrivial  solutions.  We  prove  here  that  such  irregular 
frequencies  do  not  exist. 

3.  Uniqueness 

result  can  be  stated  as  follows: 
f(0)  ♦  Ofp”1)  1  as  p  -*•  ®  , 


Specifically  our  centra 

eikop 

Theorem:  If  (a)  <P  ■  - - 


and  (b) 


♦  (p) 


♦  /  ♦(q)  ds  +  /  <Mq)  +  ky]ds 

C  *nq  ^  ®nq 


0 


for  all  p  c  CQ  and  , 
(c)  <p  is  continuous  on  CQ  u 


then  $ (p)  =  0. 


Proofs  The  proof  of  this  theorem  depends  on  the  growth  of 
potentials  with  densities  satisfying  conditions  (a) ,  (b) , 
and  (c)  of  the  theorem.  Assume  that  $  is  a  function 
satisfying  (a) ,  (b)  and  (c)  of  the  theorem  and  define  the 
functions  u+  and  u_  in  D+  and  D_  respectively  as 


'  $  (q)  (p,q)ds  +/  <f>  (q)  (p,q)  +ky  (p,q)  ]  ds  ,  {  .  ( 

u-  Co  9nq  q  Cf  3nq  °  peD_ 

As  will  be  seen  shortly,  an  essential  ingredient  involves  the 
growth  of  u_  for  large  radial  distances  from  the  origin.  Observe 
that  since  y  has  no  singularities  when  qeCQ  u  Cf  ,  p€D_  and  y 
is  a  solution  of  Laplace's  equation  it  follows  that 

V2u_  -  0  ,  p«D_.  ( 

The  jump  conditions  for  the  double  layer  defined  on  CQ  u  Cf 
take  the  form 

P*VCf  I  ♦<«>  <P'«J>d*q  •  +  *<p>+  j  *«!>  Ifc  <p,q)ds. 

pcD.  C  uC.  q  C  UC,  q  q 

t  O  X  o  r 

peCouCf  •  ( 

This,  together  with  the  continuity  of  the  single  layer,  implies 
that 

lim  3y  3y 

^ouCf  u_  (p)  «<j>  (p)  +  (q)  yjj—  (p,q)dsg+  (q)  +ky]  dsg  .  (12) 

p«D_  Co  g  Cf  q 


But  4  satisfies  the  homogeneous  equation  (b)  hence 


lim 

P+CouCf 


peO 


u  (p) 


0  . 


(13) 


However,  as  established  in  the  appendix,  lim  u_  ■  0.  Hence 

the  maximum  principle,  which  asserts  that  u_  assumes  its 
maximum  and  minimum  values  on  the  boundary,  implies  that 

u_  =  0  ,  p  e  D_  .  (14) 


Therefore 

3u-  *  0  on  C  and  C*  ,  (15) 

7rT  of 


where  indicates  the  normal  derivative  from  D_  .  Using 

the  defining  equation  (9)  for  u_  we  find  with  the  usual  jump 
conditions  for  the  single  layer 

t—  f  <Mq)  4j-(p,q)ds  +k  / <p  (q)  4--  y(p,q)ds  +  B(p)4>(p;  *  0 
P  CQuCf  dnq  q  Cf  dnp  q 


where 


B(p)  »  0  ,  p  c  Cq 

■  -k,  p  c  Cf  . 

Note  that  while  existence  of  the  normal  derivative  of  the 
double  layer  in  some  weak  sense  was  needed  in  order  to  apply 
the  divergence  theorem,  once  it  is  established  that  u_  =  0 
and  hence  has  an  ordinary  normal  derivative,  namely  zero,  the 
defining  equation  for  u_  ensures  that  the  normal  derivative 
of  the  double  layer  exists  in  the  ordinary  sense  since  u_ 
and  the  single  layer  have  ordinary  normal  derivatives. 
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Now  examine  the  limiting  values  of  u+  as  p  approaches 

C  uC-  from  D  .  Using  the  usual  jump  conditions  we  find 
o  r  + 

u+(p)*  -4(p)  +  ^P'<3^dsg+  ^q**ln""  +k7]ds  '  PeC0uCf 

c0  q  cf  * 

and,  since  $  satisfies  the  integral  equation  (b) , 
u+(p)  =  -2<J>  (p)  ,  peC0uCf  . 

Observe  that  since  4  is  assumed  to  have  growth  as  specified  in 
(aj ,  equation  18)  ensures  that  u+(p)  has  the  same  growth  on 


(17) 


(18) 


Now  form  the  normal  derivative  of  u+  from  D+  obtaining 

3-*“  *  ~>r  -  /  4  (q)  rX  (P,<3)ds  +k  /4  (q)-g^j— (p,q)  ds  -B(p)4(p)  •  (19) 

+  V  CQuCf  3nq  q  Cf  dnp  q 

Since  the  normal  derivatives  of  the  double  layer  with  continuous 
density  are  the  same  from  either  side  provided  one  of  them  exists, 
we  use  equations  (16)  and  (18)  to  obtain 

3u+ 

-28  (p)  4  (p)  =  6  (p)  u+  .  (20) 

With  the  definition  of  B(p)  (cf.  eqn.  (16))  we  see  that 
3u+ 

■  0  *  P«co'  (21) 

^  « 

afid 


3u+ 

"  "ku+  '  PeCf  * 


(22) 


[• 


since  this  property  is  inherited  from  y(p,q).  Furthermore 
by  its  construction  u+  satisfies  Laplace's  equation  in  D+ 
and  since  u+  also  satisfies  the  Neumann  condition  on  Cg  and 
the  free  surface  condition  on  C^,  u+  has  the  representation, 
following  John  (1950) , 


u .  *  E  u  (p:fl)cosh  k  (y+h)  ,  p  =  oc  +  >  a  , 

n=0  n 

where  are  the  roots  of  the  transcendental  equation  (4) 
and  a  is  any  number  greater  than  the  diameter  of  the  ship 


hull  i.e. 


a  >  max  p  . 


Recall  that  (p,9  ,y)  are  the  cylindrical  coordinates  of  the  point 
p.  Moreover,  as  shown  in  the  Appendix,  u+  «  0  ( — j^j) 


hence 


u+(P*9#h)  cosh  kn<y+h)dy  =  0  C -j^_ j) 


which  implies,  with  the  orthogonality  of  (cosh  kn(y+h)} 


on  l>2  (— h ,  0 )  ,  that 

u_  (p,3  )  =0  (-err) 


(26 


This  in  turn  implies  that 

2*  “in0 
Vp'8)  e  de 

-  0 


and  since  the  most  general  form  of  u  (p,9)  is 

ao  n 

V0-31  -  ^tannH  fi|  <V>+bn»Hfi|'  <  V>  ’  ^  < 

it  follows  that 

anm  H  fm|  <V>  »bmn  H  [m|  ,<V>  *  °  < 

Here  are  Hankel  functions  of  the  first  and  second  kind 

respectively.  The  fact  that  kR  is  positive  imaginary  for 
n>0  then  ensures  that 


Then 

u+  (p,6,o) 


Z 

n-0 


Z 


(1) 

m 


(k  p)  eimecosh  k  h 
n  n 


00 

+  *  bon  H{ml  <kob)  eim3  cosh  k  h 

m“-°°  1  1 

and  because  u+(p,  9 ,  o)  has  the  same  asymptotic  growth  as 

(hQP)  '  c^e<3n*  (18),  uniformly  in  9  we  may  conclude 

that  bQm  *  0  which  then  implies  that 
00  00 

u  (P,  9  ,y)  «  Z  Z  a  H  R (k  p)  eim  3  cosh  k  (y+h) 
n-0  m— »  ™  I  m  |  n  n 

hence  u+  satisfies  the  radiation  condition  for  -h<v<0  . 


( 


12 


Thus  u+  is  a  solution  of  the  homogeneous  floating  body 
problem  in  D+  <cf.  (1)  and  (2))  and  there fore ,  provided  that 
CQ  satisfies  the  geometric  restrictions  of  the  uniqueness 
proof  (John  (1950),  Kleinman  (1982)),  it  follows  that  u+  «  0 
in  D+  and  hence  also  on  CQuCf.  Equation  (20)  then  ensures 
that  <fr(p)  *  0  on  CQuC^.  That  is,  the  only  solution  of  the 
integral  equation  (b)  satisfying  (a)  and  (c)  is  $  »  0.  This 
means  that  the  integral  equation  (7)  has  no  irregular  frequencies 
and  has  at  most  one  solution.  Existence  of  this  solution  for 
all  k  will  be  discussed  elsewhere. 

We  remark  that  if  the  integral  equation  (7)  has  a  solution 
<J>  on  CQuCf  then  the  solution  of  the  inhomogeneous  floating  body 
problem  (1)  is  given  by 

♦  (p)  -  -i  /  *(q)  (p,q)  ds  -  \  f  ♦(q)[|j-+  ky(p,q)]ds  (3Z) 

2  3nq  q  2  Cf  3nq  q 

+  -  /  V (q)  y  (p,q)  dsq 
2  Co 

for  p«D+  u  CB> 
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Appendix :  On  the  growth  of  u+. 

Here  we  prove  the  Lemma  needed  in  establishing 
uniqueness  of  solutions  of  the  integral  equation.  For 
convenience  we  restate  it  as  follows: 

Lemma:  If  a)  $  *  —  (f  (e)  +  0(p  ) )  as  p  ♦  •  , 

/p 

b)  ♦  (p>  +  /  ♦  (q)  |“  dsq  +  /  $  (q)  [|~  +  ky]dsq  *  0  , 

PeCouCf' 

c)  ♦  is  continuous  on  CQuCf  , 


d)  u+  -  /  *(q)  ll-  dsq  +  /  ♦  (q)  tfj-  +  ky]dsa  , 


P  €  D+ 


then  u_  -  0  ( — •  as  rp  •  ,  5  <  1/2 


u+  *  0(— I72=T)  as  p] 


where  r  ■  | p | 


p  %■ y  * 

P  *P 


Proof  With  y  as  defined  in  equation  (5)  it  is  clear  that 


'  ♦(‘”!n-  dsq  ■  0  <73> 
vo  ^  p 


(A.l 


C.nB.  *<q)  +  ky]ds_  -  0(4-)  , 


(A. 2 


where  C-nB  is  that  portion  of  the  free  surface  contained 

I  a 

in  the  ball  of  radius  a.  It  is  a  bit  more  work  to  establish 


the  growth  of 


'c  ♦<*>'&"  +k^a»c 
nB,  q 


cfnBa 


where  B®  is  the  complement  of  the  ball.  Considering  first 
the  term  involving  the  normal  derivative,  which  on  Cf  is 

3  m  3 

3n  ~  by 

q  Jq 

yQ  -  0, 

we  find  that  4 


3*  - 

|  ♦  «»  a»q  -  jrj  f  ♦<«> 

nBf  0  a 


Y  +2h 


CfnBa 


R(0) 


p  dp  de 


R  (h) 


where  ( p ,  0 )  are  the  cylindrical  coordinates  of  q  on  C-  and 


R(h)  -  *r(xp-xq)‘+(zp-2q)‘+(yp+2h)‘  . 

Introduce  two  sets  of  spherical  coordinates  of  the  form 


r  sin  a  cos  e 
P  P 

r  sin  a  sin  and 
P  P 

r  cos  a 
P 


r'  sin  o'  cos  0^ 
r'  sin  a'  sin  0 


y_  +  2h  ■  r'  cos  a' 
P 


where  O<0  <2ir,  0<a<ir,  0<a’<ir/;2,  r  ,r'*0  . 

—  p—  —  —  —  —  p  — 

Clearly  (rp,0p,a)  are  the  usual  spherical  coordinates  while 
r'  and  o'  will  depend  on  h.  Explicitly 


I 
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for  p«  D+uD_.  Note  that  yp>-h  when  p«  D+uD_  hence 
o'  <*/2  whereas  o  varies  over  a  larger  interval,  in  fact 
o>ir/2  when  pcD+.  In  this  notation 


R(h)  »  /r  +p  -2r  p  sin  a  cos(0-0  ) 

P 


(A. 4) 


*<q)  dsq  *  17  |  |  ♦(q>  “ 


r  cos  a 


VB« 


0  a 


p  +p  -2rp(sin  a  cos  (e-ep)  /2 


r'  cos  a 


r'  +?^-2r*  P  sin  a'  cos  (0-0  )  3/2 


pdpdO 


(A. 5) 


It  suffices  to  consider  the  first  integral  on  the  right,  the 
analysis  for  the  second  being  identical  with  r',  a',  y'  replacing 
*p,  y*  For  brevity  we  omit  the  subscript  and  denote  rp  by  r 


in  the  ensuing  analysis  and  consider 

2ir  «■ 


r  COS  a  (  - 

I  !<r 


_ _ 4(g)  p  dp  d0 _ 

+p  -2rp  sin  a  cos  (0-9  ) ) 

P 


for  0  <  o  <  it  .  (A. 6) 


Using  the  asymptotic  form  of  4  and  the  substitution  p«rt  we  find 


2n  • 


r  cos  a 


0  a 


4(c)p  dp  do 

22  "  1 

(r  +p  -2rp  sin  a  cos(0-0_))  * 


2  n  •  ik.r t 

(  r  _  0 


cos  a  e  0  - 

77T  ^  Ci+t2- 


(f  (0)+O(l))tdtde 


2t  sin  a  cos  (0-0  )) 

P 


hence  the  integral  is  0  for  a  ?  j  .  Note  that  this  expression 

does  not  obviously  exist  when  a  -  J  .  To  see  what  happens  as  a  -  1 


observe  that 
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2ir  - 


lim  r  cos  a 
ot-*-ir/2+ 


0  a 

2ir  » 


»(q)p  dp  de _ 

(r^+p^-2rp  sin  a  cos  (6-9.J ) 


1  lim 

- 

i  j 

yp-*-0+  0  a 

1  lim  +  2ir 

ZSi  < 

(q) 


3 


*V  <(VXq)!+ Vyq)2+(VZq)2)l/2 


pdcd 


v° 


♦(q)  wr  V?'q,<lsq  ’ 


+2tr*  (p)  ,  pp>a 


where  y, 


2«<<xp-xq)2+(yp  -  yq) 2+(zp“zq) 2) 1/2 


and  the  jump-condition  for  a  double  layer  is  used.  Here  we  aake 

no  use  of  the  assumption  that  |  is  a  solution  of  the  integral 

equation  b) .  The  integral  in  the  jump  condition  vanishes  for  p 

on  Cf,  (yp-0).  Now  we  use  a)  which  asserts  that  on  Cf,  $  is 

assumed  to  grow  as  0(— ,  which  is  the  desired  growth.  Hence 

o' 

the  integral  (A.  6)  is  0  0  <_  o .  <_  it  .  Redoing  the  analysis 

with  r’ ,  a\  y’  replacing  rn ,  a,  yn  leads  to  a  similar  result. 

«,  y  >  -h.  (A. 7) 

p  *  I 

[W  +  RTRT1  odode  *  (A,8) 


Hence  we  conclude  that 


*(<J)  In -  <P>'3>dsq  *  0(T7J)  as  r 

CfhBj  9  1 


Next  we  consider 


2v 

*(q)Y(p,q)dsq  -  -  77 

.  f  ' 


Ct< 


♦  (q) 


0  a 


-  V  V*  V  V  rV  V  *■. 

•-  1.  O.sVaWVW)  u.Y 


.  v  »  ’  *"■ 


V  * V,v« 
'<  *»1 1 


wfcaJk 


Using  the  notation  previously  introduced  and  the  asymptotic  form  of 
4  we  must  treat  integrals  of  the  form 


ikQP 


0  a 


Cf  (0)+O(i)]pdpde 
“5  5 

(r  +p  ~2rp  sin  ocos(0-0p)) 


and  a  similar  integral  with  r' ,  o'  replacing  r,  a. 
The  term  involving  0 (~)  is  easily  handled  since 
2ir  -  ikQp 

e  0  ( j)  pdpde 

So  (r2+p2-2rp  sin  a  cos  (6-6  ) ) 1/^2 
0  a  P 


2  it  • 


dpde 

/pT(r  +p  -2rp  sin  o  cos  0) 


0  a 


0  0 


_ dtde _ 

St  (l+t2-2t  cos  6) 


where  c  is  independent  of  r  and  o.  This  is  seen  to  be  0 since 


the  integral  on  the  right  exists  and  is  independent  of  r. 
The  remaining  integral  is  of  the  form 
2it  ®  iknp 

_  f  f  e _ p  f(9)dpde _ 


|  j (r  +p  -2rp  sin  a  cos(9-9p)) 

0  a 


ik0a 

f  (0)  /a  e  d0 


^  (r  +a  -2ar  sin  a  cos(0-0p)) 


ik0p  d 


£+p  -2rp  sin  acos (0-9_) ) 


which  we  write  as 


I1  "  0  ^~T72^  "  — ^“1 
1  r1/2  2ikflrJ 


f  (0) 


1 

* 

.  ? 


ik.rt  2 

e  0  (1-t2 ) dt 

- - — - — 

(1+t  -2t  sin  a  cos  (0-ep)) 


2  ik-.ru  5 

+  172  %  (1'U  )dU - 3/2  >  «  • 

£  u1/2(l+u  -2u  sin  a  cos  (9-9p))J/  J 


Letting  u  *  ^  in  the  second  integral  we  get 


2ir 

I,  =  0(-fr*) - ^-77  f(9)  -T- 

1  r1/2  2ikQr1/2  |  t1' 


ik  * 

iknrt  1K0t  , 

(e  0  -e  ) (1-t  jdt _ 

2 (l+t2-2t  sin  a  cos  (0-9  )) 3/2 

r 


The  integral  in  (A. 12)  which  we  denote  as  I2  satisfies  the 
inequality 


(A. 12) 


r2  TT  1  -  W 

C  f(0)  f  (e  0 

.1 


ik„rt  ik0t 


-e  ) (1-t*) dt 
2t  sin  o  cos  (0-0p)) 


2ir  1 


7  (  *~t 


ik0rt  *kt  ^k0rt  ifeot 
-e  e  -e 


(1-t  )dtd9 


0  1/2 


2t  sin  a  cos  0) 


(A. 13) 


for  arbitrary  6  e  (0,1)  (we  further  restrict  <5  subsequently)  and 
using  the  estimates 
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<2/7,i<t<l, 
/t  ~  2 


ikflrt  ik0|  <  2 
e  — e  — 

I  , ik^rl 


ikOrt_  °t  <  4  k  r(l-t)  ,  i  <  t  <  1, 


we  obtain 


I2  1  c^lfl I.  r< 


2  it  1 


(1-t) 1+5dtde 

- 2 - - 

J (1+t  -2t  sin  a  cos  6)) 

1/2 


where  c  is  independent  of  r,  a,  ep  and  f. 

But  for  0  <_  a  <_  ir  and  0  <_  0  <_  2tt  we  may  show  that 

_ 1 _  <  _ 2  <  _2 _ 

l+t2-2t  sin  a  cos  0  l+t2-2t  cos  6  (1-t)2 

hence 


(A. 14) 


2 ir  1 


0  1/2 


dtde _ 6 

(l+t2-2t  cos  e)1”! 


The  kernel  is  weakly  singular  at  t*l,  0*0  hence  the  integral 
exists.  Thus  there  is  a  constant,  c2,  such  that 

*2  -  °2llfH-r5 

which  with  (A. 12)  establishes  that 


:x  -  0(r6  1/2)  . 


(A. 15) 


We  may  choose  <5  e  (0,  -j)  to  ensure  that  1^  decays  with  r. 
A  similar  growth  estimate  is  obtained  if  r',  o'  replace  r 
hence,  with  we  see  that 

♦ (q) Y (p,q) dSg  *  0(r5  "  1/ 2) 


This  result  taken  together  with  (A. 7)  ensures  that 


CfnBa 


4>(q)  [|~  +  kv ] ds  *  0(r6  1//2) 

q  q 


which  with  (A.l)  and  (A. 2)  establishes  that 


u+  -  0(rp  6  “  1/2) 


which  implies,  for  -h  <  u  <  0,  that 
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